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ABSTRACT 

Stellar population synthesis (SPS) models are a key ingredient of many galaxy evolution studies. 
Unfortunately, the models are still poorly calibrated for certain stellar evolution stages. Of particu- 
lar concern is the treatment of the thermally-pulsing asymptotic giant branch (TP-AGB) phase, as 
different implementations lead to systematic differences in derived galaxy properties. Post-starburst 
galaxies are a promising calibration sample, as TP-AGB stars are thought to be most prominently 
visible during this phase. Here, we use post-starburst galaxies in the NEWFIRM medium-band sur- 
vey to assess different SPS models. The available photometry allows the selection of a homogeneous 
and well-defined sample of 62 post-starburst galaxies at 0.7 < z < 2.0, from which we construct a 
well-sampled composite spectral energy distribution (SED) over the range 1200 — 40 000 A. The SED 
is well fit by the Bruzual & Chariot SPS models, while the Maraston models do not reproduce the 
rest-frame optical and near-infrared parts of the SED simultaneously. When the fitting is restricted 
to A < 6000 A, the Maraston models overpredict the near-infrared luminosity, implying that these 
models give too much weight to TP-AGB stars. Using the flexible SPS models by Conroy et al, and 
assuming solar metallicity, we find that the contribution of TP-AGB stars to the integrated SED is a 
factor of ~ 3 lower than predicted by the latest Padova TP-AGB models. Whether this is due to lower 
bolometric luminosities, shorter lifetimes, and/or heavy dust obscuration of TP-AGB stars remains 
to be addressed. Altogether, our data demand a low contribution from TP-AGB stars to the SED of 
post-starburst galaxies. 

Subject headings: galaxies: evolution — galaxies: stellar content — stars: AGB and post-AGB 



1. INTRODUCTION 

Galaxy evolution studies at all epochs strongly 
rely on stellar population synthesis (SPS) mod- 
els (e.g.. iTinslev fc Gum] H976T iLeitherer et al.1 H99l 
Bruzual fc Charlotl 120031 : iMarastonl 120051: iConrov et al.l 
2009). These models are key in deriving stellar masses, 
stellar population properties, and in some cases even red- 
shifts from spectra and broadband photometry. Thus, 
our current understanding of stellar p opulations and 
galaxy growth across cosmic time (e.g.. | Papovich et al.l 
200lHForster Schreiber et al.l l2004t iShapley et al.l l2005fc 
Kriek et al.l l2009al Ivan Dokkum et al.1 120101 7 and 



the evolution of the stellar mass density (e.g., 
Rudnick et al.l 120061: [ Stark et al.lfeOOTUMarchesini et alJ 
2009tlGonzalez et allboiOHLabbe et al.l2010D is critically 
dependent on SPS models. 
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Nonetheless, SPS models are still poorly calibrated for 
certain stellar evolution stages. In particular, the ther- 
mally pulsing asymptotic giant branch (TP-AGB) phase 
is a source of major discrepancy among different SPS 
models; different implementations of its treatment lead 
to large systematic differences in derived galaxy prop- 
erties. TP-AGB stars are cool giants (T < 4000 K) 
and primarily contribute to the near-infrared part of the 
spectrum. They are of low-to-intermediate stellar mass, 
and thus their relative contribution to the near-infrared 
spectrum is highest around 0.5-2.0 Gyr after a star- 
burst. Due to the short lifetime of AGB stars (~lMyr), 
the local globular cluster calibration samples are small, 
and the total luminosity, temperature, and carbon star 
fraction in the TP-AGB phase are not well calibrated. 
For example, t he T P-AGB phase is more prominent in 
the IM araston (2005) models than in lBruzual fe; Charlotl 
(I2003D. Therefore stell ar masses derived using the 
iBruzual fc Charlotl (|2003l ) models are generally higher 
than stellar masses derived u sing the IMarastonl (|2005|) 
SPS models (e.g..lWuvts et al.ll2007l: iMuzzin et al.ll2009t 
iKannappan fc Gawiserll2007Hvan der Wei et al.ll2006P ). 

Several studies use galaxy spectral energy distribu- 
tions (SEDs) to assess the diffe rent SPS model s , and 
resolve the discrep ancies (e.g, IMaraston et al.l 120061 : 
lEminian et al.1 120081 ). This exercise is complicated due 
to the large variety of possible stellar populations and 
the degeneracies between the star formation timescale, 
age, metallicity, and dust. However, specific galaxy pop- 
ulations can be used to constrain certain stellar evolu- 
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Fig. 1. — Synthetic filter curves of intermediate bandwidth used 
to select post-starburst galaxies. f/ m — £? m isolates the Balmer 
bre ak, while _B m — V m measu res t he slope redward s of the break. 
The lBruzual & Ch ariot ( 20Q3) and lMaraston] <2"005l) models for the 
same IMF (Salpctcr 1955) arc shown for ages of log (4/yr) fa 8.8, 9.0, 
and 9.2, and an e-folding time of log (r/yr) = 8.0. We also show an 
SED of a dusty star-forming galaxy (r = 1 Gyr and t = 100 Myr) 
with the same (/ m — -B m color as the 1 Gyr old post-starburst 
galaxy (orange). All three filters fall in the regime where the model 
colors of post-starburst galaxies are nearly identical, and thus our 
selection is not biased toward any particular SPS model. 

tion phases. For example, iConrov fc Gunnl (|2010l) use 
post-starburst galaxies in the Sloan Digital Sky Sur- 
vey (SDSS) to assess the treatment of the TP-AGB 
phase of different SPS models, as TP-AGB stars dom- 
inate the near-infrared lumi nosity during this phase . 
The y find that the mode ls bv lBruzual fc Chariot! (|2003f ) 
and IConrov et all (|2009() can reproduce the optical and 
near-infra r ed co lors of post-starburst galaxies, while the 
IMarastoiil (|2005D models cannot. 

In this Letter, we take a similar approach as 
IConrov fc Gunnl ((20101 ) and use a post-starburst galaxy 
populati on in the NEWFIRM medium-band survey 
(NMBS; Ivan Dokkum et al.ll2009l) to test different SPS 
models, and obtain new constraints on the SED shape 
during the time that the TP-AGB stars are thought to be 
most dominant. The NMBS allows us to photometrically 
select a large and clean sample of post-starburst galax- 
ies, and to construct a comp osite SED which i s mor e 
thoroughly sampled than in IConrov fc Gunnl (2010). 
Throughout the Letter, we assume a ACDM cosmology 
with Q m = 0.3, n A = 0.7, and H = 70 km s" 1 Mpc" 1 . 

2. SAMPLE SELECTION 

The NMBS uses five custom near-infrared medium- 
band filter s and covers a tota l area of 0.4 deg 2 in the 
COS MOS (jScoville et al.1120071 ) and AEGIS (IDavis et al.l 
|2007f ) fields. The medium-band filter set, in combi- 
nation with deep optical medium and broadband pho- 
tometry and IRAC imaging, provides accurate pho- 
tometric redshifts and stellar population properties 



Fig. 2. — Rest-frame C/ m — B m vs. B m — Vm for all z > 0.5 
galaxies in the NMBS (with (S/N) A - > 25). The black box indi- 
cates our post-starburst selection, characterized by red Um — B m 
colors and relatively blue Bm — V m colors. The dust vector and 
dotted line indicate that contamination by obscured star-forming 
galaxies is minimal . The galaxies are color coded following their 
Bruzual & Chariot (2003) average stellar ages. In the gray box, 
we show the Maraston (2005) average stellar ages of the selected 
galaxies. For both models the rest-frame near-infrared is down- 
w eighted when fitt i ng. T he black diamonds are the seven galaxies 
in lMaraston et al.l (120061 ). 

(e.g.. iBrammer et all 120091: Ivan Dokkum et al.1 120101 : 
IWhitaker et al.ll2010D . A full overview of this survey is 
given in K. E. Whitaker et al. (2010, in preparation). 
The photometric redshifts an d stellar population pr oper- 
ties ar e derived using EAZY ([B rammer et~al1l2008l ) and 
FAST (|Kriek et al.ll2009bf l . respectively. 

Our post-starburst galaxy sample is selected using 
three optimized synthetic rest-frame filters of interme- 
diate bandwidth^ (see Figure [1]). The U m and B m fil- 
ters isolate the Balmer break, while B m — V m measures 
the slope of the SED just redwards of the Balmer break. 
The location of the V m filter is a trade-off between be- 
ing red enough to constrain the SED slope and blue 
enough not to enter the regime where the stellar evo- 
luti on models and spectra l libraries start to deviate (see 
also iMaraston et al.ll2009() . Post-starburst galaxies have 
strong Balmer breaks, characterized by red U m — B m col- 
ors and blue B m — V m colors, and thus are expected to lie 
in the black selection box in Figure [2j The boundaries of 
the selection box are a trade-off between minimizing con- 
tamination by dusty star-forming galaxies (estimated to 
be ~1%) and obtaining a large enough sample. As indi- 
cated by the orange SED in Figure [1] dusty star- forming 
galaxies with the same U m — B m color as post-starburst 
galaxies will have redder B m — V m colors. 

In order to solely include high-quality SEDs with rest- 
frame UV-to-NIR wavelength coverage, we require a red- 

9 Rest-fra me colors are compu ted from the best-fit EAZY tem- 
plates (see Brammcr et al. 2009). 
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Fig. 3. — Composite post-starburst SED (black dots), derived from the individual SEDs (gray dots) of 62 post-starburst galaxies in the 
NMBS. The individual SEDs are normalized at 5000 A before combining. The best-fit iBTiizual fc Chariot! 120031 1 and lMarastonl 120051 1 SPS 
models to the full spectrum, with the same resolution as the data, are shown by the orange and purple curves, respectively. The red and 
blue curves show the respective best fits when excluding A > 6000 A in the fit. The inset represents the minimum \ 2 value per degree of 
freedom vs. the mean stellar age for both SPS models. In the top panel, we show the fractiona l offset between the models and the data. 
The gray shaded area indicates the uncertainties o n the composite SED. The[Bruzual & Chariot (2003]) models more accurately reproduce 
the SED shape of post-starburst galaxies than the Maraston (2005) models. 



shift of z > 0.5 and a if-band signal-to- noise ratio (S/N) 
of 25. We also exclude galaxies for which the IRAC fluxes 
are strongly contaminated by nearby sources, an d thus 
cannot properly be deblended (jLabbe et al J 12005). The 
final selection consists of 62 galaxies at 0.695 < z < 
2.028, with a median redshift of z — 1.56. Six targets 
have spectroscopic redshifts, yielding a median photo- 
metric redshift uncertainty of Az/(1 + z) = 0.005. 

We have inspected the SEDs of all targets by eye and 
none have obvious problems. Furthermore, for none of 
the selected galaxies there are indications that the rest- 
frame UV-to-NIR continuum is dominated by dusty star 
formation (Figure EJ. Similar to spectroscopic post- 
starburst galaxy samp les, obscured star formation ma y 
have been missed (e.g.- lMiller fc Owenl[200ltlGotoll2007f ). 
However, some contamination by a dusty star-forming 
population will not compromise our results (see Section 
4.1). 



To obtain a homogeneous sample and to avoid mod- 
eling degeneracies, our selection intentionally excludes 
younger or heavily dust-obscured post-starburst galaxies 
and specifically targets galaxies with short star forma- 
tion timescales, for which the TP-AGB phase is most 
prominent. Therefore, our sample is likely more restric- 
tive than other spec troscopically selected post-starburst 
galaxy samples (e.g.. lQuintero et al~ll2004h . However, for 
the purpose of this Letter it is not relevant whether our 
sample is complete, as long as the selection is not biased 
toward any particular SPS model. 

3. COMPOSITE SPECTRUM 

In order to make a composite spectrum, we deredshift 
the SEDs of all galaxies to rest frame, without changing 
the observed fluxes. All SEDs are normalized at 5000 A, 
the longest wavelength for which the variations between 
the different SPS models are negligible (see Figure [1]) . 
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TABLE 1 

Stellar. Population Properties 







logt 


logr 


A v 




log Af« 


log SFR log SFR/M, 


ID 


Fit region Model 


(yr) 


(yr) 


(mag) 


Z 


Mq 


^©yT -1 yr" 1 X? cd 



Composite SED Full range 



Ma-3650 

Zspec = 1.91 



A < 6000 A 



Full range 



BC03 

M05 

BC03 

M05 

BC03 

M05 



9.06 
9.05 
9.06 
9.06 
9.08 
8.94 



o.oi 

0.01 
0.00 
0.00 
+0.18 
0.02 
+0.01 
0.01 
+0.07 
0.19 
+0.11 
0.17 



8.10 
8 
8.10 
8.10 
8.10 



+0.04 
0.02 

10 +ooo 



7.90 



0.00 
+0.20 
0.05 
+ 0.08 
0.03 
+0.14 
1.10 
+0.22 
0.90 



0.16 
0.00 
0.16 
0.08 
0.40 
0.10 



+0.02 
0.07 
+0.00 
0.00 
+0.15 
0.12 
+0.08 
0.08 
+0.45 
0.13 
+0.13 
0.10 



0.020 

0.020 
0.020 
0.020 
0.008 
0.020 



11 io+ - 00 

±1.±U_Q Q2 

in QQ+0- 00 
iu.y»_ 00 

11 in+ 009 

11.1U_q.05 

11 q+ 003 
i±.uy_ 12 

10.761°;?? 



-0.63 
-0.67 
-0.63 
-0.65 
-1.03 
-1.54 



+0.01 
0.10 
+0.00 
0.00 
+0.18 
0.13 
+0.10 
0.12 
+0.47 

27.34 
+0.73 
20.98 



7 o+0.00 
LL - ,o -0.08 

-11 66+ 00 
i±.oo_ 00 

11 7S+ - 14 
1A - 1 J -0.08 

-11 74+ - 20 

1J -- '^-0.08 
11 Q9+ - 42 
+0.70 
-20.97 



-12.30 



0.73 
5.08 
0.72 
0.94 
2.91 
2.90 



Note. — We fit a grid with log(i/yr) between 7.6 and 10.1, in steps of 0.01 (not exceeding the age of the universe), log(r/yr) between 7 
and 10 in steps of 0.1, and Ay between and 3 in steps of 0.02 mag. For each model we explore three different mctallicitics: subsolar (0.01 
and .008 f or the [M araston (2005) and Bruzual & Chariot (2003) models, respectively), solar, or supersolar (0.04 and 0.05 for the IMarastonl 
]2005l) and IBruzual & C hariot (200__ models, respectively). To account for systematic errors, we increase all photometric uncertainties by 
quadratically adding 5% of the photometric flux. The errors represent the 68% confidence intervals. 



We create a composite SED by averaging the rest-frame 
fluxes in wavelength bins of 20 data points each. The 
uncertainties on the composite spectrum are determined 
by bootstrapping. 

Figure [3] shows the composite post-starburst galaxy 
SED in black. The photometric data points of the in- 
dividual galaxies are represented by the gray dots. The 
uniformity of all individual SEDs is remarkable and re- 
sults in a high-quality, low-resolution spectrum of a post- 
starburst galaxy. Even subtle features, such as the Mg n 
absorption line at 2800 A and the continuum break at 
2640 A, are detected. 

4. COMPARISON TO MODELS 

4.1. Wruzual & Charloi 1(200 3 ) and \Marastwl K200A) 

The high quality of the composite post-starburst 
galaxy SED allows the assessment of the different SPS 
models. We fit the c omp osite SEP by b oth the 
IBruzual fc Charlotl (|2003l) and IMarastonl (|2005l ) models, 
assuming a star formation history (SFH) of the form 
ip{t) tx exp(— f/r) and leaving age (t), the e-folding 
time (r), the amount of dust attenuation (Ay), and 
metallicity (Z) as free parame t ers (s ee note to Table [T]). 
We assume the ICalzetti et al.l ((2000h attenuation curve, 
which we imple ment as a uniform screen. We adopt the 
iSalpeterl (|1955f ) initial mass function (IMF), as t his IMF 
is avai labl e for both SPS m odels. Compared to a lKroupal 
(pfffll or IChabrierl (pOOl IMF our choice will primar- 
ily affect the mass-to-light ratio and has little impact 
on other stellar population properties. For each bin, we 
determine the effective filter curve, by adding the dered- 
shifted filter curves of all included data points. In most 
cases, this is a variety of broadband and medium-band 
filters with different observed wavelengths, as our sample 
spans a large range in redshift. 

We s tart by fitting the fu ll wave lengt h range. T h e best 
fits for IBruzual fc Charlotl pOOl and IMarastonl (2005) 
are represented by the orange and purple fits i n Fig- 
ure [3] respectively. The IBruzual fc C hariot (2003) mod- 
els yie ld an acceptable fit to the data, while the lMarastonl 
(2005) models do not fit the full wavelength range simul- 
taneously. The best-fit stellar population properties are 
broadly consistent, with a slig htly lower value for stel- 
lar mass (M*) and Ay for the IMarastonl (|2005l) models 
(Table [J). We have explored other SFHs as well, such 
as a delayed exponential SFH (ip(t) ~ t exp(— t/r)), 



a truncated SFH preceded by a constant star forma- 
tion rate (SFR, with r = t st0 p — i s tart), a truncated 
SFH preceded by an exponentially increasing SFR, and 
a two-component population, but none prov ide a sig- 
nifican tly better correspondence between the IMarastonl 
r2()()5; models and the composite SED. The SFH used to 
generate the model SEDs in Figure [3] produces a sharply 
peaked spectrum; adding an older or dusty star-forming 
population would only broaden the model SED, which 
wo uld then produce an even larger discrepancy between 
the [Maraston] (|2005h models and the data. 

To better understand the discrepancies between the 
models in the rest-frame near-infrared, we repeat the fit- 
ting, but now restricting the wavelength range to A < 
6000 A. The IBruzual fe Charlotl (|200l and IMarastonl 
(2005) models, as represented by the red and blue curves, 
respectively, fit this wavelength range nearly equally well 
and yield consistent v alues for M*, t, r, and A y (Ta- 
ble [_). However, while IBruzual fc Charlotl (120031) also fit 
the full wavelength coverage, the IMarastonl ([2005) mod- 
els overpredict the rest-frame near-infrared flux. 

4.2. Tests and Comparison to other Studies 

As a check on our procedure, we stack the best fits to 
the individual SEDs for both models. The stacks are sim- 
ilar to the purple and orange curves in Figure __ Along 
the same lines, we assess th e qualities of the individ - 
ual fits. All galaxies favor the IBruzual fc Charlotl (|2003l ) 

models, with median values for Xbco3 an< ^ Xmo5 0i? 0.91 
and 4.3, respectively. 

Wherea s bot h the IBruzual fc CharlotJ (|2003l) and 
Maraston ( 2005) SPS models can reproduce the B m — V m 
and C/ m — B m colors in the selection box, the exact col- 
ors at fixed age and star formation timescale are slightly 
different. In order to test whether these small differences 
cause a bias that favors a particular SPS model, we split 
the post-starburst galaxy sample at B m — V m = 0.55, and 
repeated the analysis for both sub-s amples. Both SEDs 
are si gnificantly better fit by the IBruzual fc Charlotl 
(2003) models (with similar fit qualities as for the full 
composite spectrum), where the redder sample is older 
by ^0.1 dex. Additionally, we repeat the analysis for 
a redder boundary B m — V m < 0.6, yielding a sample 
of ^110 galaxies. The composite spectrum is better fit 
by a slight ly older stellar population , but still strongly 
favors the IBruzual fc Charlotl ([2003) models above the 
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Fig. 4. — Comparison of the composite post-starburst galaxy spectrum with the FSPS models by Conroy ct al. (2009), for which the 
temperature and bolometric luminosity of the TP-AGB stars can be modified with respect to the Padova stellar evolution models (Alog T 
and AlogL, respectively). The left panel illustrates the influence of varying the luminosity, while keeping the temperature fixed to the 
best-fit value. In the right panel, the luminosity is fixed to the best-fit value, but the temperature varies. In both cases, the best fit is 
shown in purple. For both panels, t, t, and Ay are fixed as well. The inset shows the best-fit AlogT and AlogL (purple square), and its 
corresponding 68% and 95% confidence intervals. 



IMarastoril (|2005f h Altogether, this confirms our unbiased 
s election. 

iMacArthur eTail (poloh compare the near-infrared 
fluxes of two spiral galaxies to predictions fro m long-slit 
spectr a, reporting better co nsistency with the | Marastonl 
(2005) models than with IBruzual fc Chariot! (|2003j. 
However, the interpretation of their results is compli- 
cated by the presence of a significantly older stellar pop- 
ulation, which contributes 26%-67% to the total light, 
and the presence of dust, which may mimic the effects of 
a TP-A GB star populati on. Nonetheless, it is possible 
that the IMarastoril (|2005h model indeed provides a bet- 
ter description for the SED shape during the evolution 
p hase of these spi r al galaxies. 

iMaraston et all pOM 12001 use the broadband pho- 
tometry of seven spectroscopically confirmed galaxies to 
assess the fit quality of different SPS models and find no 
overall preference for a particular SPS model. Only one 
galaxy (ID 3650) of this sample would have entered our 
selection (see Figure [2]). We have re- fitted this galaxy 
using the same method as for the composite spectru m, 
and - broadly consistent with IMaraston et al.l ((2006) 
we neither find a strong preference for a particular SPS 
model (see Tabled]). Thus, the differences in fitting tech- 
niques do not account for the large discrepancy in % 2 
value for the composite spectrum. Moreover, as several 
galaxies in our sample, in partic ular at higher re dshift, 
are almost equally well fit by the IMaraston! (I2005D mod- 
els, ou r study is not in disagreement with lMaraston et aLl 
((20061) . 

This redshift dependence may be caused by the fact 
that the metallicity and d ust content at fixed stellar mass 
evolv e with redshift (e.g-. lErb et al.l l2006: Maiolin o et al.1 
2008), and thus it may be possible that the contribution 



from TP-AGB stars also changes with time. We test this 
by splitting our sample at z — 1.56 and repeating the 
analyse s for both sub-samp l es. A lthough, for both sam- 
ples the lBruzual fc Char iot (2003) models pro vide signif- 
icantly better fits, th e relative fit quality of thelMarastonl 
(2005) models versus lBruzual fc Chariot! (|2003| ) increases 
with redshift. Future work might be able to separate 
these trends and constrain TP-AGB models as a func- 
tion of stellar mass and metallicity. 

4.3. Utilizing a Flexible SPS model 

The difference in fit quality betwe en the 
iBruzual fc Chariot! (|2003l) and IMaraston! (|2005l ) models 
for the composite post-starburst galaxy spectrum reflects 
the different treatments of the TP -AGB phase. Our 
result s suggest that the treatment bv IBruzual fc Chariot! 
( 2003) is more appropria te for our post-sta rburst galaxy 
sample than that of the Maraston (2005[) models. We 
use the flexible SP9ZB (FSPS) models by IConrov et al.l 
(200^ 120101) to obtain a more quantitative constraint on 
the TP-AGB phase, as FSPS allows the modification of 
the bolometric luminosity and effective temperature of 
the TP-AGB stars. 

We start with the best-fit FSPS model to the com- 
posite post-starburst galaxy spectrum, restricting the 
wavelength region to A < 6000 A and assuming the 
latest default Padova TP-AGB models dGirardi et al.1 
[2lM[Marigo fc Girardil[2007l: IMarigo et al.ll2008[ ). Simi- 
lar to the other models, this is a stellar population with 
log (i/yr)=9.04, log (r/yr)=8.1, and an Ay of 0.4 mag 
(for solar metallicity). Next, we fix the age, star forma- 

10 These models can be downloaded at www. cfa. harvard. edu/~ 
cconroy/FSPS.html. 
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tion timescale and dust content of the stellar population 
and vary both the effective temperature and bolomet- 
ric luminosity of the TP-AGB stars, quantified as shifts 
with respect to the Padova evolutionary tracks, AlogT 
and Alog L, respectively. Thus, we ignore any potential 
degeneracies between AlogT, AlogL, and other stellar 
population properties (jConrov et al.ll2009l ). 

Figure [4] illustrates the influence of both parame- 
ters on the SED and shows their x 2 contours. The 
best fit has a reduced y 2 val ue of 0.75, thus simi- 
lar to iBruzual fe Chariot! ((2003). The composite post- 
starburst galaxy spectrum suggests that the predicted 
effective temperatures of TP-AGB stars in the Padova 
models are consistent with our observations, but the 
overall luminosity is ^0.5 dex lower. 

5. CONCLUSIONS 

In this Letter, we define a photometrically selected 
sample of 62 post-starburst galaxies from the NMBS and 
use their composite SED to obtain new constraints on 
the SED shape during the time that the TP-AGB stars 
are th ought to be most d o minan t. The SED is well fit 
bv the IBruzual fc Charlotl (|2003i ) SPS models, while the 
iMarastonl (|2005f T models do not reproduce the rest-frame 
optical and near-infrared parts of the SED simultane- 
ously, implying that these models give too much weight 
to TP-AGB stars. Th is has previously been found by 
iConrov fc Gun n (2010) using post-starburst galaxies in 
the SDSS. 

The high-resolution photometric sampling of the 
NMBS allows us to derive quantitative constraints on 
the lumino sity in the TP-AGB phase. Using the FSPS 
models bv IConrov et al.l ([2009') . we find that the bolo- 
metric luminosity of TP-AGB stars is a factor of ~3 
lower than predicted by the latest Padova TP-AGB mod- 
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els, when assuming solar metallicity. However, an inde- 
pendent abundance measurement is needed to break the 
degeneracy between metallicity and the TP-AGB phase 
parameters. 

The significant reduction in the bolometric luminos- 
ity of TP-AGB stars that is required to fit the observed 
post-starburst SED reflects one or more failures of the 
current generation of SPS models. This reduction can 
be physically achieved in the models by reducing TP- 
AGB lifetimes, reducing bolometric luminosities of stars 
in the TP-AGB phase, and/or embedding a significant 
fraction of TP-AGB stars within optically thick circum- 
stellar dust shells. These different explanations may 
be disentangled by appealing to mid-IR data (see, e.g., 
IKelson fc HcJde^l20ia ISalim et al.l[2009h . 

Taking our results at face value, we infer that stellar 
masses an d other st e llar p opulation parameters derived 
using the Maraston] (|2005[ ) might be biased for certain 
evolutionary phases. We stress that our findings do not 
imply that the Bru zual fc Charlotl (|2003l ) models are best 
for all galaxies or for all evolutionary phases. Using simi- 
lar techniques as presented in this Letter, in combination 
with flexible SPS models, which allow unrestricted mod- 
ification of uncertain evolution phases and input param- 
eters (such as metallicity) , it will be possible to substan- 
tially reduce many of the major uncertainties plaguing 
all SPS models. 
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